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ABSTRACT 

In this paper we refine the calculation of primordial helium recombination, accounting for sev- 
eral additional effects that were neglected or treated more approximately in previous studies. 
These include consideration of (i) time-dependent radiative transfer interaction between the 
2'Pi - 1 1 So and 2 3 P; — 1 1 So resonances; (ii) time-dependent radiative transfer for the partially 
overlapping n'Pi - 1 1 So, n l T>2 — 1 So and n 3 Pi - l'So series with 3 < n < 10; (iii) electron 
scattering within a kernel approach. We also briefly discuss the effect of electron scattering 
and H i quadrupole lines on the recombination of hydrogen. Although the physics of all con- 
sidered processes is interesting and subtle, for the standard cosmology, with helium abundance 
F p 0.24, the overall correction to the ionization history during helium recombination with 
respect to the previous implementation of CosmoRec remains smaller than \AN e /N e \ - 0.05%. 
The dominant improvement is caused by consistent inclusion of resonance scattering for the 
2*Pi - 1 1 So- For cosmologies with a large helium fraction, Y p ^ 0.4, the difference reaches 
lAA^e/A^I — 0.22% at z - 1800, however, the overall correction to the CMB power spectra is 
small, exceeding \ACi/Ci\ - 0.05% only at I > 3000. In comparison to the current version of 
Recfast the difference reaches |AQ/Q] - 0.4% at / =* 3000 for Y p 0.4, and also for the 
standard value Y p 0.24 we find differences \ ACi/Ci\ > 0.1% at / 2500. The new processes 
are now included by the cosmological recombination code CosmoRec and can be activated as 
needed for most settings without affecting its runtime significantly. 

Key words: Cosmic Microwave Background: cosmological recombination, temperature 
anisotropies, radiative transfer 



1 INTRODUCTION 



Over the past fe w years, several independent groups have recon- 
sidere d ( e.g. see iDubrovich & Grachev! 120051 : Chluba & Sunyaev 



2006): 
2008al: 



Kholupenko & IvanchikI 12006): Hirata 

Wong&Scottl 120071; iRubino-Martm et alj [2 008: 
Karshenboim & Ivanovl 120081 : iHiratal 120081 : (Chluba & Sunyaev 



20081 : Ijentschural 120091; ILabzowskv et al J 120091 ; IGrin & Hir~ 
20ld : lAli-Haimoud & Hirata! 1 20 id) the cosmological recombina- 
tion problem of helium and hydrogen, with the premise of assessing 
how uncertainties in the theoretical treatment of the recombination 
process could affect the science return of Planck and future CMB 
experiments. Neglecting these r ecently-derived cor r ection s to the 
standard recombination model l lSeageretalJll999l. [2 000 ) biases 
cosmological parame ter estimates teubino-Martm et"al] l20ld : 
IShaw & Chliibal 1201 lh , in particular, the inference of the slope 
of primordial scalar perturbations and baryon content. Especially 
when trying to distinguish between different competing inflation- 
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ary models (see iKomatsu et al"]|201 lL for recent constraints from 
WMAP), such inaccuracies should be avoided. 

These efforts are now coming to an end, and it appears that all 
important corrections to t he standard recombination scenario have 
been identified (e.g. see iFendTe t al. 2009; Rubino-Mart in et al] 
[20Kl for overview). This has lead to the devel opment of two new, 
indep endent reco mbination codes, CosmoRec JChluba & ThomasI 
1201 ll) and HyRec I Ali-Hai'moud & Hiratal fcOllh . both of which in- 
clude all important modifications to the recombination problem, 
superseding the physical model of Recfast. These new codes are 
fast, accurate, and self-contained in the sense that they do not 
need externally-calibrated effective rate modulation ("fudge") fac- 
tors from more complete recombination codes. A detailed code- 
comparison is currently in preparation, however, there are a few 
subtle radiative transfer effects during helium recombination that 
have not been discussed in detail thus far. 

The current helium recombination mod el of CosmoRec is 
mainly based on the approach presented in IChluba & Sunyaev! 
< T2010d) . which we briefly recap in Sect. 12.21 Here we extend this 
model by explicitly solving the radiative transfer problem for the 
evolution of the high frequency distortion of the CMB introduced 
by neutral helium resonances. This extension allows us to include 
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redistribution of photons by resonance and electron scattering, and 
the full time-dependence of photon production and absorption (see 
Sect, l2,lt in an efficient and accurate manner. 

The refinements to the helium recombination model discussed 
here are especially interesting since current constraints on the he- 
lium mas s fraction, Y B , deriv ed using ACT l lDunklev etaT]|201lh 
and SPT l lKeisler et alj|201 lh data indicate a slightly larger helium 
abundance of Y p = 0.313 ± 0.044 and Y p = 0.296 ± 0.03, respec- 
tively. In that case it is expected that, e.g., the effect of resonance 
scattering on the 2 1 P ] - 1 1 So line is modified as compared to the 
standard cosmological model, with Y p = 0.24. As this effect is cur- 
rently only taken into account approximately (i.e., by using a cor- 
rection functio n for the escape probability which was calibrated at 
fixed Y p = 0.24 lRubino-Martm et al.l2008h . one might worry about 
the precision of the helium recombination calculation. Further, a 
detailed radiative transfer computation treats the feedback and in- 
teraction of overlapping excited levels (see Sect. 12. 2t . Our calcu- 
lations show that for the analysis of future data from Planck, the 
considered processes are only marginally important (see Sect. |3.5l l. 
Nevertheless, the new version of CosmoRec (version 2.0) now in- 
cludes a helium radiative transfer module that can be activated as 
neededQ This module might be useful for recalibration of the Rec- 
fast fudge-factors (see Sect. I3.5.21 . once the final recombination 
code comparison is completed. 

Finally, we briefly consider the effect of H i quadrupole lines 
during hydrogen recombination. So far, only the rcd-ls transitions 
were inc luded using an effect ive transition rate between the nd and 
np state dGrin & Hiratal201fj) . Here we confirm this approximation 
using a full radiative transfer computation. Furthermore, we include 
the quadrupole transitions among the excited levels and show that 
these actually dominate the correction. However, the final effect is 
negligible (see Sect. l3.7t . 



2 HELIUM RECOMBINATION PHYSICS 

In this section we briefly outline how the radiative transfer calcula- 
tions during helium recombination are carried out. The approach 
strong ly relies on the treatment presented in IChluba & Thomas] 
d201 lh for the case of hydrogen recombination with several im- 
portant differences, as we explain here. 



2.1 Time-dependent radiative transfer: main aspects 

In simple words the helium and hydrogen recombination problem 
can be described as follows: free electrons are captured into some 
excitecQ. atomic state, forming a neutral atom under the emission 
of a photon in the free-bound continuum. From the excited lev- 
els, electrons can cascade towards energetically lower states, re- 
leasing one or more resonant photons. While most of the possible 
sequences of transitions are compensated by their inverse process 
(mediated by interactions with photons drawn from the ambient 
CMB), as recombination proceeds, more electrons reach the ground 
state of the atoms and distortions to the CMB blackbody spectrum 
are created. These distortions are very small in the Rayleigh- Jeans 
tail and close to the maximum of the CMB blackbody, however, 
in the Wien-tail they exceed the thermal radiation field (e.g. see 



1 CosmoRec is available at www . Chluba . de/CosmoRec 

2 Direct recombinations to the ground state can be neglected 
IZeldovich et alJl968l ; |Peeblesll968tlctiluba & Sunvaevll2007l) 



Table 1. The first few main transitions to the ground-state of helium. The 
resonances are ordered according to their transition frequency, v,i s , cor- 
responding to the sequence of radiative feedback between lines. We give 
the Az/z — Av/v necessary to reach the next lower-lying resonance. For 
comparison, due to the thermal motions of the atoms one has Avd /v =: 
1.6 X 10~ 5 [z/2000] 1 ' 2 . The line center resonance cross section at redshift 



z — 2000 is given by rx r 



1.3 X 10 24 w, A,i s v.. 3 , with w, 



for the n Pi - 1 1 So and n 3 Pi - 1 1 So series, and vv, = 5 for the singlet 
quadrupole lines. The last column also indicates the distance from the line 
center in Doppler units up to which the lines are optically thick at z — 2000. 



Level 


A,i s [s- 1 ] 


Viis [PHz] 


cr r [cm 2 


1 


|Az/z| 




thick 
D 


2 3 P! 


177.58 


5.0691 


5.1 x 10" 


21 






0.63 


2'Pi 


1.799 x 10 9 


5.1305 


5.0 x 10" 


14 


1.2 x 10 


-2 


120 


3 3 P! 


56.1 


5.5631 


1.2 x 10" 


21 


7.8 x 10 


-2 




3 ! D 2 


1.298 x 10 3 


5.5793 


4.7 x 10" 


20 


2.9 x 10 


-3 


2 


3'P! 


5.663 x 10 8 


5.5824 


1.2 x 10" 


14 


5.6 x 10 


-4 


33 


4 3 P! 


23.7 


5.7325 


4.8 x 10" 


22 


2.6 x 10 


-2 




4'D 2 


748.5 


5.7394 


2.5 x 10" 


20 


1.2 x 10 


-3 


1.4 


4 l Pi 


2.436 x 10 8 


5.7408 


4.9 x 10' 


15 


2.4 x 10 


-4 


13 


5 3 P, 


12.1 


5.8100 


2.3 x 10" 


22 


1.2 x 10 


-2 




5'D 2 


431.4 


5.8135 


1.4 x 10" 


20 


6.0 x ur 4 - 


1.2 


5'P, 


1.258 x 10 8 


5.8143 


2.4 x 10" 


15 


1.4 x 10 


-4 


7 



ISunvaev & Chlubal2009l. for more details about spectral distortions 
from the recombination epoch). In particular, non-thermal photons 
emitted at high frequencies (for hydrogen by the Lyman- series and 
for helium the n'Pj - 1 'So, n 3 Pi — 1 1 So and «'D2- l'So resonances) 
can lead to re-excitation of electrons from the ground state to higher 
levels and the continuum, and hence inhibit the recombination pro- 
cess. Therefore the recombination problem requires simultaneously 
solving the evolution of the populations of electrons in the excited 
states as well as the high frequency spectral distortions introduced 
by the recombination process. In particular it is important to com- 
pute at which rate photons can leave the cores of the resonances, 
thereby escaping from further (strong) interactions with the atoms. 

Without additional approximations the recombination prob- 
lem described above is prohibitively complicated. However, in 
the case of hydrogen it has now become possible to treat the re- 
combination process with sufficient detail, achieving a precision 
< 0.1% in the ionization history of the Universe close t o z— 1 100 
dChluba & Thomas! 1201 ll; lAli-Haimoud & Hiradl201 ll) . Here we 
extend the approach o flChluba& Thomas! feOllh to account for the 
detailed physics of the helium recombination process. The main ad- 
dition is that we explicitly solve the radiative transfer problem for 
the He i resonances, instead of relying on simple (occasionally even 
rough) approximations. 

For the line emission and absor ption process we use the 1 + 1 
photon description (see Sect. 3.2 in lchluba&Thomasll201ll) . We 
neglect two-photon corrections to the shapes of the line profiles, as 
these are expected to lead to small changes in the ionization history 
dSwitzer & Hiratal2008bf) . Nevertheless, the 1 + 1 photon treatment 
allows accounting for the full time-dependent evolution of the pho- 
ton distribution, respecting detailed balance in the distant wings 
of the line profiles. Both aspects have been sh own to be very im- 
portan t for hy drogen recombination (e.g., see IChluba & Sunvaevl 
l2009bl . l2010bl) . but thus far were only partially treated for helium. 

We also include the redistribution of photons over frequency. 
Close to the Doppler cores of the optically thick lines, resonance 
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scattering is most important. Here we implement a Fokker-Planck 
approximation for the redistribution functior^, however, in addition 
we compare the results w ith those from a Monte Carlo approach 
dSwitzer & Hirat dl2008 j) . In the distant wings of the resonances 
and also in the cores of the optically thin lines, redistribution of 
photons by electron scattering becomes significant. We explain the 
treatment of this process in more detail below (see Sect. 12. 3b . 

Finally, for helium recombination there is one very im- 
portant difference relative to the hydrogen recombination prob- 
lem: since already at z — 2000 a small fraction of neutral 
hydrogen atoms is present, high frequency photons emitted in 
the main helium resonances are able to additionally ionize hy- 
drogen atoms. This removes helium photons from the cores 
of the lines and consequently strong l y accelerates hel i um re- 
combination dKholupenko et alj |2007| ; Iswitzer & Hiratal |2008eJ ; 
iRubifio-Martfn et alj 120081) . Using the variables x = V/V21 and 
An x = v 2l An v , where v 2l = 5.1305 x 10 15 Hz is the 2 x Pi - 1% 
transition frequency, and A/i v denotes the deviation of the photon 
occupation number from the CMB blackbody, the corresponding 
term in the transfer equation reads 



dAn x 



dz 



of t '(x)N?,'c 



An x , 



(1) 



Here cr^'(x) is the photoionization cross section from the ground 
state of hydrogen and iVj? 1 denotes the Hi ground state popula- 
tion. This expression has to be added to the radiative transfer equa- 
tion describing resonance scattering, and line emission/absorption 
for the individual helium resonances (see Table [T] for transition 
ra tes and frequencies for so me of the neutral helium lines) given 
in lChluba & Thomas! (EoTTh . 

In the work presented here we do not explicitly discuss correc- 
tions to the 2'So - 1 1 So two-photon channel. For example, stimu- 
lated 2 1 So- 1 1 So two-photon decay should be included, accelerating 
helium recombination; at the same time one expects that feedback 
of photons from the main helium resonances can re-excite atoms 
to the 2' So by 1 1 So — 2' So two-photon absorption. Similar pro- 
cesses a re very important for the 2s- Is two-photon channel of hy - 
drogen dChluba & Sunvaevll2006l ; iKholupenko & Ivanchikl 12009) . 
but for He i recombination the 1 1 So - 2 'So two-photon channel is 
much less impo rtant, allowing only ^ 8% of the helium atoms to 
become neutral jRubino-Martm etal]|2008l) . We implemented the 
corresponding problem, and found that a detailed treatment leads 
to a correction < 0.01% in the ioni zation history of helium, w hich 
is in agreement with the analysis of lSwitzer & Hirat3 d2008bl) . 



2.2 Overlapping resonances in the upper shells 

In the previous treatments of He i recombination, the feedback of 
photons from higher to lower lines is modelled by the approx- 
imation that all the resonances are isolated (see Fig. [TJ, in the 
sense that within the proximity of any individual line, only the 
line itself contributes sign ificant opacity dSwitzer & Hirata 2008al ; 
IChluba & Sunvaevll2010al) . In this approach the standard Sobolev 
escape probability, Psi^) = [1 - e^s]/-^, where is the Sobolev 
optical depth of the i' h resonance, is replaced by the effective escape 



3 The terms for the radiative transfer equation are given by Eq. (9) and (11) 
in IChluba & T homas (2011), with appropriate replacements related to the 
differences between hydrogen and neutral helium. 



Case 1 




Isolated 
line 



Feedback photons 




Case 2 




Overlapping 
lines 



Feedback photons 




b c 



Figure 1. This figure illustrates the concepts of "isolated" and "overlap- 
ping" resonances. In Case 1 the resonances "a" and "b" are isolated, and 
the only interaction between them is mediated by redshifted photons from 
"b" to "a". In Case 2, the resonances "a" and "b/c" can be treated as isolated 
lines, while "b" and "c" feel the presence of each other. 



probability (see lChluba & Sunva ev 2010a, for more details) 



pi ~ p< 



Anl _ r c 
1 e i- 



Art 



P' 



1 - »' P' 

1 r™ ' d 



(2a) 
(2b) 



Here, P\ = Ps(p' d T' s ), where p' d = 1 - p' m is the death proba- 
bility or line albedo of the i th resonance; An' L (z) is the deviation 
of the line occupation number of the i* resonance from a black- 
body; An'_(z) e J-* is the distortion of the photon occupation num- 
ber caused by the energetically higher line j at the line center of i 
after redshifting and absorption by neutral hydrogen; finally, . 
is the Hi continuum opacity between the upper and lower line. 

In the case of hydrogen (TjL^ = 0), the isolated line approx- 
imation works very well for the Lyman-series up to n - 10 - 
20, since the resonances are sufficiently separated in frequency, 
while for larger n, overlap of the opticall y thick resonances be- 
comes important |Ali-Hai'moud et alj|2010h . As pointed out earlier 
dChluba & SunvaevteOloX for the resonances in the upper shells 
of helium, the isolated line approximation should already break 
down for n > 3. This is because in every shell with n > 3 pho- 
tons encounter a triplet of n 3 P t — 1 1 So, rc'D 2 - 1 1 So and n'Pj - l'So 
resonances (see Fig.[T](. In particular, the zi'D2 - 1 1 So quadrupole 
and n'Pi - 1 'So dipole lines are separated by < 30 Doppler widths 
at z — 2000, while the n [ ~Pi - l'So dipole lines are optically thick 
up to comparable distances from the line center (see Table Q] for 
examples). Such proximity between the lines allows for a cross- 
communication in the damping wings, which becomes increasingly 
important towards the end of helium recombination; therefore, to 
correctly evaluate the escape probability of each line, the photon 
fields of nearby resonances must be taken into account. Further- 
more, redistribution of photons by electron and line scattering may 
lead to additional subtle changes in the recombination dynamics 
that are not captured by the simple isolated line approximation. Ex- 
plicitly solving the full radiative transfer equation, as described in 
the previous section, accounts for the physics of this problem self- 
consistently. We explain the different effects in more detail below. 
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2.3 Electron scattering 

To account for the effect of electron scattering on the cosmologi- 
cal recombination process, we start with the Boltzmann collision 
term describing the redistri bution of photons over freq uency in an 
isotropic mediurrQ (e.g., see lSazonov & Sunvaevll200(J) : 

^ = rP^ATjlAns^-^-AnJdV. (3) 

<5T e e _ sc Jo 1 

Here An v is the deviation of the photon occupation number from 
the CMB blackbody; t c = jo-jN^cdt is the optical depth with 
respect to Thomson scattering; cr x ss 6.65 x 10~ 25 cm 2 is the Thom- 
son cross sections; N c is the free electron number density; and 
Vy^y'iTz) denotes the electron scattering kernel for thermal elec- 
tron at temperature T e . For the computations presented here, we 
us e the Kompaneets equation kernel, Pk, as defined by Eq. (23) 
ofi azonov & Sunvaevl d2000h . It accounts for Doppler broaden- 
ing, Doppler boosting, and electron recoil. However, both electron 
recoil and Doppler boosting remain rather small, and the redistri- 
bution of photons by electron scattering is dominated by Doppler 
broadening. 

When the variations in the solution for the spectral distortion 
are very broad in comparison to the characteristic width of the scat- 
tering kernel, it is possible to resort to a Fokker-Planck approxima- 
tion of the Boltzmann collision integral, Eq. l[3}. In this approach, 
Any is approximated using the first and second derivative of the 
photon occupation number, and the integral over frequency can be 
carried out analytically, leading to a second order parabol ic par- 
tial differential equation (e.g., see lSazonov & Sunvaevll2000h . Such 
types of problems can be efficiently solved numerically, and the 
cosmological recombination code, CosmoRec, is readily extended 
to take this into account. 

However, as explained earlier l lHirata &~Sw itzer 2008]), the 
characteristic width of the electron scattering kernel during helium 
recombination is given by 

£„J^ MWx i<r*JS«J^£. (4) 

v V m c c 2 V 2000 V m c V 

where Avq c is the Doppler width of the helium resonances, and 
V«?He/'«e ~ 86 arises from the helium to electron mass ratio. This 
turns out to be much wider than the typical variation of the spec- 
tral distortions found on the blue side of the high frequency Hei 
resonances, and hence the simple Fokker-Planck approximation is 
expected to break down in regions close to the resonances. In this 
case, one faces solving the numerically expensive dense matrix 
equation M u = b for the evolution of the photon occupation num- 
ber during helium recombination. 

Fortunately, it turns out that in places where the variation of 
the phase space density is large (i.e. in the vicinity of the Doppler 
cores of the main Hei resonances), electron scattering is always 
much less likely than any of the line processes. This can be seen by 
comparing the mean free path to Thomson scattering, Z T = 1 /W c <x T , 
with the one for resonance interactions (scattering + absorption) 
at the line center, /, ^ l/Nn E1 cr r . Here Ahei is the neutral helium 
number density and cr, the line center resonance scattering cross 
sections, which at z — 2000 is given in Table Q] for the first few 
lines. At z - 2000 one has N He[ /N h =s 0.04, so that even for the 
5 3 Pi - 1 1 So resonance (the one with the longest mean free path 



among the lines given in TableQJ we find Ij/h a 14. For the singlet 
n'Pi - 1 1 So resonances one finds li/l t » 1 even at very early stages 
of neutral helium recombination. Similarly, lj/l r > 1 for the shown 
quadrupole lines, with a comfortable margin. 

One can therefore resort to an iterative approach when solv- 
ing the dense matrix equation at fixed redshift z: given the matrix 
equation Mo u = b describing the solution for the photon occupa- 
tion number without accounting for the effect of electron scatter- 
ing, one can compute the solution u , benefiting from the extreme 
sparseness of the problem. This solution is then used to compute 
the correction vector Sb = —SM uq, which can be derived from the 
Boltzmann integral for electron scattering^ In the next iteration, 
one now has to compute the correction to the solution, Suo, from 
the matrix equation M du = db , which again can be calculated 
efficiently using banded matrix routines. This procedure can be re- 
peated until convergence is reached, however, we find that usually 
one iteration per time-step is sufficient. 

The abo ve procedure is prac t ically equivalent to the approach 
explained by lAli-Ha'imoud et all d2010h in the case of hydrogen 
recombination. However, during hydrogen recombination electron 
scattering is much less important for three reasons: (i) the width 
of the electron scattering kernel is only about 86 VWT^Hc - 43 
times larger than the Doppler widths of the H i resonances; (ii) all 
the H i Lyman-series lines are extremely optically thick, so that in 
the vicinity of the r esonances line scattering al ways dominates over 
electron scattering dChluba & Sunvaevll20093) : (iii) the number of 
free electrons rapidly decreases as hydrogen recombination pro- 
ceeds. Overall, electron scatteri ng never becomes very im portant 
during hydrogen recombination dChluba & Sunva ev 2009a). How- 
ever, even in the case of hydrogen recombination, the effect of elec- 
tron scattering on the recombination dy namics requires a treatmen t 
using the precise redistribution kernel Jab -Haimoud et al. l l2O10l) . 
We briefly discuss this in Sect. 13.31 

We emphasize that in most astrophysical situations the redis- 
tribution of photons by electron scattering can b e safely described 
using a Fokker-Planck approximation (e.g., see IPozdniakov et"al] 
1 19791) . The cosmological recombination epoch is one example for 
which this approximation fails, as we show below. It is the inter- 
play between line emission/absorption processes, causing narrow 
features in the photon distribution, and broad electron scattering 
which render it important to use the full redistribution kernel. 



3 MODIFICATIONS TO THE COSMOLOGICAL 
IONIZATION HISTORY 

In this section, we describe the corrections to the cosmological re- 
combination history of helium arising from the additional processes 
treated in this paper. We proceed in a step-by-step manner, first 
comparing the results of our treatment with those obtained in ear- 
lier works. As a fiducial model we use the cosmological parameters 
fi b = 0.044, Qdm = 0.216, T = 2.725 K, H = 71.0kms-* Mpc" 1 , 
Y p = 0.24, and N v = 3.046. In addition, we assume flatness and an 
equation of state parameter w = - 1 for dark energy. 



4 We neglect stimulated electron scattering and only follow the spectral 
distortion with respect to the CMB blackbody (see Chluba & Thomas 20 1 1 , 
for further justification). 



5 In essence, one has to insert the solution no into Eq. {3) and carry out the 
integrals to determine bo- Since the integral has to be computed for every 
grid point, this computation still remains rather time-consuming. 
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Figure 2. Effective escape probability of the 2 3 Pj — 1 1 So resonance for dif- 
ferent settings, regarding feedback. For comparison, we show the standard 
Sobolev escape probability. 

3.1 Transport in the 2' Pi - 1 ' S and 2 3 Pi - 1 1 So resonances 

In this section, we consider physical processes in the 2 'Pi - 1 1 So 
and 2 3 Pi - 1 1 So transitions which most strongly regulate the forma- 
tion of the neutral helium. One way to represent the corrections to 
the recombination dynamics of each line is by computin g the effec- 
tiv e escape probabili t ies (seelChluba & Sunvaevl ( feOlOal) . Sect. 4.1; 
or lSwitzer&Hiratal fe008ah . Sect. V). This can be achieved in a 
multi-step process: (i) calculate the solution to the recombination 
problem based on an approximate analytic prescription for the es- 
cape probabilities of each line; (ii) compute the escape probability 
with a radiative transfer code using the obtained solutions for the 
level populations; (iii) iterate until convergence is achieved. Start- 
ing with the prescription implemented in CosmoRec v 1.3 we found 
that usually one iteration with the new helium radiative transfer 
module was sufficient. Below we present the results for several 
cases, comparing with the standard Sobolev approximation and dif- 
ferent intermediate approximations. 
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2'p, - l's Singlet line 




Figure 3. Effective escape probability of the 2 Pi - 1 So resonance. The 
2 1 P i - 1 1 So resonance is treated as isolated line, i.e. no feedback is included. 

the 2'Pi - 1 1 So resonance) no longer reach the 2 3 Pi - l'So in- 
tercombination line because of absorption in the Hi Lyman con- 
tinuum. This feedback model was already incorporated to Cos- 
moRec vl.3, however, the escape probability is slightly smaller 
when using the results from the radiative transfer code of Cos- 
moRec v2.0 (curve labeled 'radiative transfer (with feedback)'). At 
early times differences at the level of ^ 10% are visible, while at 
z 2000 the effective escape probability again becomes identical 
to the case without feedback^ This result for the escape proba- 
bility is consistent wit h the simpler radiative trans fer computation 
presented in Fig. 6 of IChluba & Sunvaevl d2010ah . Physically the 
main addition of the new radiative transfer module discussed here 
is line scattering, which can be neglected for 2 3 P t — l'So- However, 
the feedback model of CosmoRec vl.3 in addition neglected part 
of the time-dependence, s ince the more detailed computation of 
IChluba & Sunvaevl fcOlOal) was too time-consuming to be solved 
quickly. CosmoRec v2.0 now includes all these effects simultane- 
ously without compromising the performance too much. 



3.1.1 Escape from the 2 3 P[ - l'So resonance 

Figure[2]presents the effective escape probability of the 2 3 P 1 — 1 1 S 
resonance for different cases. For comparison, we show the escape 
probability in the standard Sobolev approximation (red/dashed 
line). For the curve labeled 'RM2008', w e use the integral approx- 
imation from lRubino-Martih et al.l J2008I) . Appendix B, Eq. (B.3), 
to account for the increase in the escape probability by absorption 
of He i photons in the Hi Lyman continuum. Here the 2 3 Pi - 1 1 So 
resonance is treated as isolated line. The curve labeled 'radiative 
transfer (no feedback)' is obtained with the new helium radia- 
tive transfer code of CosmoRec v2.0, again assuming that only 
the 2 3 Pi - l'So resonance is present. The differences relative to 
'RM2008' are rather small, and only visible close to z — 1900, in- 
dicating that with the radiative transfer treatment escape from the 
2 3 Pi - 1 So resonance is slightly slower. 

If we now include the feedback among all the levels up to 
n = 10, we obtain the other two curves shown in Fig. [2] The curve 
la beled 'CS2010' is computed using the simplified feedback model 
of IChluba & Sunvaevl ( [2010au . in which all helium resonances are 
treated as independent distant lines. Feedback stops being impor- 
tant at z < 2200, since photons from the upper lines (mainly 



3.1.2 Escape from the 2 Pi - l'So resonance 

Figure [3] shows the isolated line computations for the 2'Pi - l'S 
resonance. We ran two cases with the radiative transfer code of 
CosmoRec v2.0, one with (curve labeled 'radiative transfer (no 
feedback)') and another without line scattering ('radiative transfer 
(no feedback, no diffusion)'). Again we compare with the standard 
Sobolev approximation, as well as the escape probability used in 
Cosmo Rec vl.3 (curve labeled 'RM 2008'), which is based on the 
work of lRubifio-Martfn et alj d2008l) . 

At low redshifts, the effect of H i continuum absorption be- 
comes very large, leading to a strong increase i n the escape prob- 
ability of the 2' Pi — l'Sp line (c ompare with ISwitzer & Hiratal 
l2008al ; iRubino-Martm et"alll2008h . Like for the 2 3 Pi - l'So in- 
tercombination line, this is accounted for using a ID integral ap - 
proximation (Appendix B, Eq. (B.3) of lRubino-Martm et al.l2008l) . 
but in addition with a redshift-dependent correction function to 
account for the effect of line-diffusion on the escape probability. 

6 To compute the escape probabilities in all cases, we use the solution for 
the populations of the levels that did not include the feedback correction. 
This allows us to directly compare all of them in one figure. 
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Figure 4. Modification to the effective escape probability of the 2' P[ - 1 1 So 
resonance caused by interaction with the other high frequency lines. 



feedback from the upper resonances with n > 2. Accounting for 
those lines, we find changes at the level of a few percent with re- 
spect to the isolated line case (see curve labeled 'interaction with 
all lines n < 11'). Obviously, here the main effect arises from the 
feedback of photons from the triplet of lines with n = 3. Again one 
can observe that at z < 2300, feedback ceases since photons are 
efficiently absorbed in the H i Lyman continuum before being able 
to reach the 2'Pi - l'So line. 

In Fig. [4] for comparison we also show the modification in the 
escape probability obtained with the feedback model of CosmoRec 
vl.3, w hich is based on the approximations of IChluba & Sunvaevl 
d2010al) . Although qualitatively similar, the feedback term is under- 
estimated in this model. This is likely caused by an overestimation 
of the Hi continuum opacity between the resonances of the third 
and second shell, as in CosmoRec vl.3 it is only computed roughly, 
using the Kramers' approximation for the cross section. Also, the 
red wing absorption of 2 3 P i - 1 1 So photons by the 2 1 P i - 1 1 So line 
is not accounted for. With the new radiative transfer module these 
approximation can be avoided. 



In lRubifio-Martm et alj J2008h this effect was calibrated for fixed 
Y p = 0.24, using a simple radiative transfer calculation, while here 
we consistently treat the transport problem case-by-case. As Fig. [3] 
indicates, for Y v = 0.24 at both high an d low redshifts, good agree- 
ment with iRubino-Martfn et all d2008h is found. The main differ- 
ences are visible close to z — 2300, reaching the level of 1% - 5%. 
At this redshift the optical depth of the resonance is largest, so 
that similar to the escape problem of Hi Lyman-o' photons (see 
IChluba & Sunvaevll2009bl. for more detailed explanation), the oc- 
cupation number in the distant wings becomes important, introduc- 
ing time-dependent aspe cts which are not included by the much 
simpler approximation of lRubino-Martfn et alj J2008t> . 

From Fig. [3] one can also clearly see that at low redshifts 
line scattering leads to an important additional increase of the es- 
cape probability. We confirmed that the result of the new radiative 
transfer model is in a greem ent with the previous computations of 
IRubino-Martfn et alj d2008T) when neglecting line-scattering. The 
effect of line-diffusion is currently not included by HyRec, but does 
lead to a correction of A/V e /7V e -0.3% at z — 1800, however, for 
the analysis of CMB data fr om Planck with Y„ ^ 0.24 this effect 
can apparently be neglected dShaw & Chlubalr201lh . 

As next step, we also include the effects of the other lines on 
the effective escape probability of the 2'P[ - 1 1 So line. If we only 
add the 2 3 Pi - 1 1 So intercombination line in the radiative transfer 
code, we obtain the solid black line in Fig. [4] For clarity, we plot 
the relative difference with respect to the case for which only the 
2'Pi - l'So resonance is present. Although the 2 3 Pi - l'So line 
energetically lies below the singlet linefj the escape probability is 
affected at the level of a -0.35% at z — 2300, a correction that 
is not accounted for in earlier feedback models. This correction is 
consistent with absorption of 2 3 Pi - l'So photons in the red wing 
of the 2' Pi - 1 1 So resonance. Since the Sobolev optical depth of the 
2 3 Pi - 1 1 So resonance during helium recombination remains Ts < 3, 
reabsorption of 2' Pi - 1 1 So photons close to the 2' Pi - l'So line 
center in the distant blue wing of the 2 3 P[ - 1 ' So line - an effect that 
could increase the 2 Pi - l'So escape probability - is negligible. 

However, the modification to the escape probability caused by 
the presence of 2 3 P[ - l'So photons is smaller than the effect of 

7 This implies that there is no direct line feedback to the 2' Pi - 1 ' So reso- 
nance by redshif'ting. 



3.2 Partially overlapping triplets of lines with n> 2 

In this section we illustrate the modifications to the escape prob- 
abilities of the upper, partially overlapping triplets of resonances 
with n > 2. Figure [5] shows the effective escape probabilities of 
the different lines for 3 < n < 5. For the intercombination lines, 
one can clearly see that th e simple isolated line approximation of 
IChluba & Sunvaevl d2010ah is in rather good agreement with the re- 
sult of the full radiative transfer calculation carried out here. Only 
at low redshifts, where the effect of H i continuum absorption is be- 
coming significant, the difference is large. Th is is because in the 
simplifi ed version of the treatment described in lChluba & Sunvaevl 
d2010ij) which is used in CosmoRec vl.3 this correction was ne- 
glected to accelerate the computation. Similarly, for the rc'Pi - 1 1 S 
singlet lines at e arly times we find very good agreement with the ap- 
proximations of IChluba & Sunvaevl d2010ah . with differences only 
showing up at the end of helium recombination. 

For the «'D2 - 1 1 So quadrupole lines the overall differences 
relative to the radiative transfer calculation are much larger. We at- 
tribute these differences to details in the radiative transfer between 
the partially overlappi ng resonances, as well as the validity of ap- 
proximations made in IChluba & Sunvaevl d2010ah . For example, 
for n > 3 the quadrupole and singlet lines are sufficiently close 
to each other to start communicating directly. Photons close to the 
rc'Pi - l'So resonance can be absorbed in the blue wing of the 
n'D 2 - l'So line. Similarly, absorption of photons in the cores of 
n'D2 - l'So lines can occur via the red wing of the n'Pi - l'So. 
The latter process is more important, since the optical depth of the 
n'Pi - l'So lines is orders of magnitude larger than for the cor- 
responding quadrupole lines. The line therefore remains optically 
thick up to larger distances from the line center, comparable to the 
distance between the lines (see TableQJ. The isolated line approxi- 
mation breaks down for these lines, and at larger n even the partial 
overlap with the intercombination line should start mattering. 

We also confirmed numerically that line scattering itself is 
only responsible for a small part of the difference. Given that the 
differences are rather small, implying small changes in the recom- 
bination history, we did not attempt to further isolate the physical 
reasons for the deviations. Overall, for the upper levels the main 
improvement of the new radiative transfer module with respect to 
CosmoRec vl.3 is the consistent inclusion of line scattering and the 
absorption of photon in the H i continuum, which both become im- 
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Figure 5. Comparison of escape probabilities for the n 3 P[ — 1 ' So inter- 
combination lines (upper panel), the n'D2 - 1 1 So quadrupole lines (middle 
panel), and the ra'Pj - 1 1 So singlet lines (lower panel) with 3 < n < 5, 
The red/dashed lines give the Sobolev approximation; t he solid/black lines 
are obt ained using the isolated line approximation of IChluba & Sunvaevl 
l2010al) . The blue/dashed-dotted lines are computed with the radiative 
transfer code of CosmoRec v2.0. 



Figure 6. Effect of electron scattering on the effective escape probability of 
the 2 3 Pi - 1 1 So intercombination line. For the first three curves the 2 3 Pi - 
1 1 So resonance is treated as isolated line, while for the latter three curves 
all resonances with n < 11 are included. 



portant only at the end of helium recombination. However, because 
the upper levels contribute much less to the total recombination 
rate, the impact of these modifications remains small. 



3.3 Electron scattering 

In this section we investigate the effect of electron scattering on the 
escape probabilities. In particular, we demonstrate that the simple 
Fokker-Planck approximation for the redistribution of photons over 
frequency is a rather poor description of the problem. Unless stated 
otherwise, for the results presented here we use the helium radiative 
transfer code of CosmoRec v2.0. 



3.3.1 Effect of electron scattering on the 2 3 Pi — 1 1 So line 

We first consider the isolated 2 3 P! - l'So intercombination reso- 
nance. Figure|6]presents the effective escape probability for differ- 
ent runs with electron scattering enabled. The solid/blue line repre- 
sents the result from Fig.[2]for comparison. The dotted/red curve is 
obtained using the Fokker-Planck approximation for the redistribu- 
tion function of electron scattering. Comparing to the case without 
electron scattering, at z ^ 2100 electron scattering causes a de- 
crease in the escape probability, while at low redshifts it allows 
more photons to escape. If we compute the effective escape proba- 
bility using the electron scattering kernel, as explained in Sect. 12.31 
we obtain the black/dot-dash-dashed curve in Fig. [6] Here one can 
see that electron scattering only affects the escape probability at 
high redshifts and that the effect is significantly smaller than in 
the Fokker-Planck approximation. This curve is also in ve r y goo d 
agreement with the result of Fig. 4 in Iffirata & Switzerl i200A 
which was computed with a Monte Carlo approach (see Sect. 13.41 
for more details), which confirms the validity of the approximations 
used in our implementation of this process. 

To better understand the cause for these differences, in Fig. [7] 
we present the spectral distortion introduced by the 2 3 Pi — 1 1 So 
intercombination resonance at z — 2395 as an example. Without 
any electron scattering, the spectrum has a steep step close to the 
resonance frequency of the 2 3 P| - 1 1 So line at x = 0.988. As the line 
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Figure 7. Effect of electron scattering on the spectral distortion (in arbitrary 
units) caused by the 2 3 Pj - I 1 So intercombination line. The inset shows a 
small region around the transition frequency. 



itself is only mildly optically thick, this curvepractically resembles 
the case of no frequency redistribution at alllfj 

If the effect of electron scattering is included using the Fokker- 
Planck treatment, we obtain the red/dashed curve. Two things are 
evident: (i) the step in the spectral distortion is strongly smeared 
out as a result of electron scattering, and (ii) the total intensity 'step' 
from the far blue to the far red side is reduced. However, it turns out 
that inside the Doppler cor^l of the resonance (see inset of Fig.[7J 
the total intensity actually increased by ^ 20%, explaining the ob- 
served decrease of the escape probability in Fig. [6] relative to the 
case without electron scattering. The main reason for this increase 
is that photons from the red side of the resonance scatter back into 
the line center a nd to the blue side of the re sonance (this process 
was discussed in IChluba & Sunvaevl j2009ah for the H i Lyman a 
resonance), thereby enhancing the line intensity in the core. This 
inhibits photon escape from the resonance, in agreement with what 
is seen in Fig. [6] 

If we now compute the solution for the 2 3 Pi — 1 1 So spectral 
distortion using the Kompaneets kernel, P K (see Sect. 12. 3| for ex- 
planations), we again observe a reduction of the total intensity step 
from the far blue to the far red side (Fig. [7] blue/dash-dotted line). 
However, this time the step is practically not smeared out. Instead 
we can see the effect of distant jumps of photons in the tail of the 
blue wing distortion. These are caused by the rare electron scat- 
tering events, when accounting for the non-diffusive nature of the 
redistribution process. Photons are kicked out of the Doppler core 
and redistributed over a wide range of frequencies instead of be- 
ing smeared out over the region close to the Doppler core, where 
photon production and absorption are most important. Also, one 
can observe a gradual reduction of the distortion in the red wing 
of the resonance, indicating that even at 100 - 200 Doppler width 



We usually call this case the "no-scattering" approximation. 
9 For the 2 3 Pi — 1 1 So intercombination line most of the support of the 
excited level by photons is coming from within a few Doppler widths of the 
resonance. This is because the 2 3 P[ - 1 1 So line is not very optically thick, so 
that absorption in the damping wing can be neglected. In contrast to this, for 
the 2' Pi - 1 1 So line even photons at separation s of a ^ 100 - 1000 Dopple r 
widths matter (for more detailed explanation see IChluba & Sunvaevl2010j) . 



below the resonance photons are efficiently scattering towards the 
blue wing, where a photon deficit is found. 

Again taking a close look at the distortion across the Doppler 
core of the line shows that the support of the 2 3 P t level by photons 
is increased with respect to the no-scattering solution, however, not 
as much as in the Fokker-Planck approximation. This explains why 
the effective escape probability is decreased by a smaller amount 
than for the Fokker-Planck case (cf. Fig.©. 

At low redshifts (z < 2000), we find that the kernel approach 
for the spectral distortion in the cases with and without electron 
scattering practically coincide in the Doppler core, while in the 
Fokker-Planck approximation we find a rather large reduction of 
the line intensity. This does explain the different behaviour in the 
escape probability (cf. Fig.|6](, however, what is the physical reason 
for this? When using the Fokker-Planck approximation, the distor- 
tion again is strongly smeared out in the vicinity of the Doppler 
core, and photons from within the Doppler core diffuse into the 
wings. This time also absorption of photons by neutral hydrogen 
becomes important, so that there are far fewer photons on the red 
side of the resonance, which in principle could return to the line 
center by electron scattering. Before they can do so they are likely 
absorbed in the H i Lyman continuum. Similarly, photons that scat- 
ter from the line center to the blue side of the resonance are likely to 
be absorbed by hydrogen before they return to the core to feed back. 
Therefore diffusion of photons out of the Doppler core wins and the 
total intensity is reduced with respect to the case without electron 
scattering. In the kernel approach, the smoothing of the spectral 
distortion is basically negligible, and the equilibrium of emission 
and absorption terms in the core remains largely unaltered by rare 
electron scattering events. 

If we finally also include the effect of the other lines in the 
computation of the effective escape probability of the 2 3 Pj - 1 1 S 
resonance, we obtain the other three curves shown in Fig. [6] Again 
we give the result without electron scattering (see Fig.|2j for com- 
parison. In the Fokker-Planck approximation at high redshifts, the 
total escape probability is slightly smaller than for the case without 
electron scattering, and at z < 2000, where the 2 3 Pi — 1 1 So res- 
onance practically behaves as an isolated line because of H i con- 
tinuum absorption, we again see a significant increase in the escape 
probability^ The overall behaviour of the escape probability in the 
Fokker-Planck approximation is in good agreement with our under- 
standing of the problem without the other lines included: at early 
times electron scattering delays recombination in comparison to the 
case without electron scattering. 

When using the kernel approach, something more surprising 
happens. Instead of additionally decreasing the escape probability 
at high redshifts, we can observe a small increase at z > 2500. 
Comparing with the previous behaviour for which no additional 
lines were present, this is not expected. One reason for this effect 
is that photons approaching the 2 3 P t — 1 1 So resonance from the 
blue side occasionally scatter directly to the red side of the Doppler 
core, essentially jumping over the line without feeding back. In the 
Fokker-Planck approximation, this process is much less important, 
since large jumps in frequency are not accounted for. However, the 
more crucial effect is actually a reduction of the escape probability 
from the 2'P[ - 1 'So singlet line once electron scattering is included 
(see Sect. |3.3.2l ; this reduces the spectral distortion on the blue side 
of the 2 3 Pi - 1 So resonance (cf. Fig. [8j( and hence decreases the 



The curves with and without feedback but electron scatter treated using 
a Fokker-Planck and kernel approach are indistinguishable at z < 2100. 
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Figure 8. Effect of electron scattering on the spectral distortion (in arbitrary 
units) in the vicinity of the 2 3 Pi - 1 1 So and 2'Pj - 1 1 So resonances. 
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Figure 9. Effect of electron scattering on the effective escape probability of 
the 2'P] - l'So singlet line. We show the relative difference in the escape 
probability with respect to the isolated line case. The blue dashed line is 
shown for comparison and only includes the feedback correction discussed 
earlier (see Fig.|4). 

feedback correction. As shown in Fig. [9] at z — 2900 the escape 
probability of the 2'Pi - 1 1 So line is reduced by ^ 14% when ac- 
counting for electron scattering in the kernel approximation. This 
implies that the feedback correction should also be reduced by a 
similar amount. We find this being in agreement with the result pre- 
sented in Fig. [6] To confirm this statement we in addition ran the 
radiative transfer code but including electron scattering only in the 
vicinity of the 2 3 Pi - 1 1 So line. In this case the feedback correction 
was much closer to the case without electron scattering. 

3.3.2 Effect of electron scattering on the 2'Pi — 1 'So singlet line 

We carried out a similar exercise for the 2' P] -l'So singlet line. The 
results are shown in Fig. [9] Let us first focus on the curve obtained 
using the Fokker-Planck approximation (red/dash-dot-dotted line). 
We can see that with respect to the case without electron scattering 
the escape probability is increased at all redshifts. At low redshifts, 
where H i continuum absorption is important, this is expected like 
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Figure 10. Effect of electron scattering on the escape probability of the 
;) 3 Pi - 1 1 So intercombination lines. 

for the 2 3 Pi - l'So intercombination line, however, escape from the 
2'Pi - 1 1 So singlet line is not affected by as much (in the case of 
the intercombination line the relative difference was ~ 3 - 4 times 
larger). Here it is important to mention that the 2'Pi - l'So sin- 
glet line itself is very optically thick, so that the spectral distortion 
is controlled by the line over a much broader range of frequen- 
cies (roughly a few hundred Doppler widths as opposed to only a 
few Doppler widths for the intercombination line). Therefore the 
effect of electron scattering on the redistribution of photons is lim- 
ited much more to the line wings. 

To understand the increase in the escape probability at high 
redshift, we ran the isolated line case for the 2'Pi - l'So resonance, 
but including electron scattering using the Fokker-Planck approach. 
This showed that the main increase in the escape probability is ac- 
tually due to modifications in the transfer problem of the line itself, 
rather than being a consequence of changes in the blue wing feed- 
back correction. Again diffusion of photons out of the line center 
into the wings is helping to reduce the total support of the 2'Pi 
level, and hence increase the escape rate. 

If we now turn to the results obtained within the kernel ap- 
proach, we can see that at high redshifts photon escape is hindered, 
while at low redshifts again an increase can be seen, but not as 
large as for the Fokker-Planck approximation. Again we ran the 
case with all lines included (solid/black line) and when treating the 
2'Pi - l'So resonance alone (red/dot-dash-dashed line). Compar- 
ing the two shows that changes in the solution for the line itself are 
causing the largest difference. At high redshifts, the main reason for 
the difference is that more photons from the red wing are scattered 
back to the Doppler core of the line and also into the blue wing 
when encountering a free electron. This delays photon escape. At 
later times, when hydrogen absorption starts becoming important, 
the extra broadening introduced by electron scattering mainly re- 
moves photons from the Doppler core, reducing the line intensity. 

3.3.3 Effect of electron scattering for the upper lines 

To close the discussion about the effect of electron scattering on the 
escape probabilities from the helium resonances, Fig.[l0]shows the 
results of our computations for the « 3 Pi - l'So intercombination 
lines with n = 3, 4 and 5. At high redshifts, like in the case of the 
2 3 P! - l'So line, we can see that the feedback correction is slightly 
smaller than in the case without electron scattering. As mentioned 
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Figure 11. Effective escape probability of the 2 Pi - 1 1 So resonance. Com- 
parison between CosmoRec v2.0 and the Monte Carlo treatment of the iso- 
lated line with and without electron scattering. 



Figure 13. Effective escape probability for the 3 1 D2 - 1 1 So quadrupole line. 
Comparison between CosmoRec v2.0 and the Monte Carlo treatment of the 
isolated line for the effect of electron scattering. 
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Figure 12. Effective escape probability of the 2 Pj - 1 So resonance. Com- 
parison between CosmoRec v2.0 and the Monte Carlo treatment of the iso- 
lated line for the effect of electron scattering. The optical depth for coherent 
line scattering is ~ 300 that of true (incoherent) absorption, so the Monte 
Carlo must draw many scattering events per trajectory, making it slower to 
integrate than the triplet or quadrupole lines (4 CPU-days per point). 

above, this is mainly because the blue wing distortion of the n 3 Pi - 
1 1 So lines is reduced by electron scattering, with photons close to 
the optically thick n}V\ — 1 So resonances returning to the Doppler 
core, where photon emission and absorption are effective. 

For the quadrupole lines we find a similar effect, while for the 
higher n'Pi - l'So singlet lines, the correction is in the opposite di- 
rection, again indicating that modifications to the radiative transfer 
problem of the line itself dominated. This is in agreement with the 
result found for the 2'Pi - 1 1 So line (see Fig. [9}. 

3.4 Comparison with Monte Carlo computation 

To confirm the results obtained with the He i radiative transfer mod- 
ule of CosmoRec, we compare the escape probabilities for several 
cases with those of a Monte Carlo implementation of the escape 
problem. The Monte Carlo follows an ensemble of photon trajec- 
tories subject to redshift, atomic and electron scattering, line ab- 



sorption, and photon destruction in the H i continuum . The me thod 
used here is described in detail in ISwitzer & Hiratal j2008al) and 
ISwitzer & Hiratd ( t2008bh . where it is shown that the probability 
distribution of photons in the Monte Carlo solves the continuous ra- 
diative transport equations. The simulation of line scattering draws 
from t he probabi l ity dis tribution for the velocity of the scattering 
atoms and finds the outgoing photon frequency 

through kinematics that assume dipole angular redistribution and 
recoil. Electron scattering is treated identically, except that the elec- 
tron velocity is normally distributed. True absorption and emission 
are simulated through the Voigt profile (with no memory of the in- 
coming photon frequency). 

Figures QT| [12] [13] compare the radiative transport and Monte 
Carlo solutions for the escape probability in the 2 3 Pi - l'So, 
2'Pi - 1 1 So, and 3'D2 - 1 'So transitions, respectively. The ordinary 
Fokker-Planck approximation breaks down for electron scattering 
because that scattering produces jumps in photon frequency which 
are both large and rare compared to scattering through the line. Be- 
cause the Monte Carlo fully simulates these scattering processes, 
it demonstrates that electron scattering is treated accurately in the 
kernel approach developed here. Agreement in transport through 
2'Pi - 1 1 So also shows that the Fokker-Planck approximation is ac- 
curate for coherent atomic scattering. The Monte Carlo is computa- 
tionally expensive, and cannot currently form the basis for a practi- 
cal recombination code, but corroborates the fast kernel method in 
CosmoRec v2.0. 

We also computed the effective escape probability for the iso- 
lated 3 'Pj - l'So singlet line finding excellent agreement with the 
Monte Carlo result. This further confirms the validity of the approx- 
imations used in the new radiative transfer module of CosmoRec 



3.5 Changes in the ionization history and power spectra 

In this section we discuss the effect of different processes on the 
dynamics of helium recombination. We start by switching off line 
and electron scattering, directly comparing CosmoRec v2.0 with 
CosmoRec vl.3. The results for the modifications in the ionization 
history for Y p = 0.24 are presented in Fig. [14] The improvement 
over CosmoRec vl.3 is mainly due to consistent inclusion of time- 
dependent aspects, line-scattering, a precise modelling of Hi con- 
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Figure 14. Direct comparison of CosmoRec v2.0 with CosmoRec vl.3. The 
recombination history of CosmoRec vl.3 was used as reference and both 
line and electron scattering are switched off. The differences are because of 
time-dependent aspects and additional approximations made in the simpli- 
fied feedback treatment of CosmoRec vl.3. 



Figure 16. Convergence of the electron scattering correction with number 
of included shells n. We internally compare CosmoRec v2.0 with and with- 
out electron scattering enabled, modelling it with the kernel approximation. 
Electron scattering leads to a deceleration of helium recombination at high 
redshifts, while at low redshifts it accelerates recombination. 
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Figure 15. Convergence of the line scattering correction with number of 
included shells n. We internally compare CosmoRec v2.0 with and without 
resonance scattering enabled. Line scattering leads to a slight deceleration 
of helium recombination at high redshifts, while at low redshifts it acceler- 
ates recombination. We also show the direct comparison of CosmoRec v2.0 
with CosmoRec vl.3 (used as reference). 

tinuum absorption, as well as the overlap of lines for the upper 
states. The differences are small and should affect the predictions 
for the CMB power spectra at a negligible level. 

In Fig.[l5]we show the separate effect of resonance scattering 
on the ionization history. Resonance scattering leads to an acceler- 
ation of helium recombination around z ~ 1800. The convergence 
with number of shells is very fast, with the most important correc- 
tion coming from n = 2. 

We can also compare directly with CosmoRec vl.3, but ac- 
counting for the effect of line scattering on t he 2'Pi - 1 1 Sp reso- 
nance using a calibrated correction function l lRubino-Martm et"al] 



|2008|) . As Fig.[T5]shows, the agreement is excellent, with improve- 
ments at the level of \AN C /N C \ - 0.04%. Such a difference leads to 
negligible changes in the CMB temperature and polarization power 
spectra and even for future CMB experiments can be omitted. How- 
ever, as we discuss below when varying the helium abundance, 
larger differences appear. 

Finally, in Fig. [Tfj] we show the effect of electron scattering 
on the helium recombination history. Again the improvement over 
CosmoRec vl.3 for this case is negligibly small, and convergence 
with n is very rapid. Given that one full run with all high frequency 
helium lines up to n = 10 included and electron scattering modelled 
in the kernel approximation takes about a 5 hours on eight CPUs, 
this is fortunate. Neglecting electron scattering, a helium atom with 
10 shells takes about 90 seconds on one CPU. Also, when electron 
scattering is modelled with the full kernel, but the correction is only 
included in the vicinity of the 2'Pi — 1 1 So and 2 3 Pi — 1 1 So lines, the 
run finishes in roughly 15-30 minutes on eight CPUs, depending 
on the total number of lines included. For numerical computations 
this simplification is very useful and we checked that the difference 
in the ionization history is negligible. The results presented below 
are calculated using this approximation. 

To conclude, for the standard cosmological model, with he- 
lium abundance Y p = 0.24, the improvement of CosmoRec v2.0 
over CosmoRec vl.3 remains smaller than \AN e /N e \ — 0.05%, and 
the difference in the CMB power spectra computed with these two 
versions of the recombination code are expected to be negligible. 

3.5.1 Dependence of the correction on the helium abundance 

As we just saw for the standard cosmology, the overall improve- 
ment of the helium recombination calculation are not very signifi- 
cant, and smaller or comparable to o ther corrections (e.g., bec ause 
of two-photon and Raman processes l lSwitzer & Hiratal2008bh dur- 
ing helium recombination), that are omitted here. The currently 
used approximations of CosmoRec vl.3 therefore already provide a 
precise description for the problem. However, when going slightly 
away from the standard cosmology the differences start increasing 
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Figure 17. Total correction to the ionization history for the standard cosmo- 
logical model with varying helium abundance. We compare the result ob- 
tained with the helium radiative transfer module of CosmoRec v2.0 with the 
one computed with the more approximate treatment of the earlier version, 
CosmoRec vl.3. The main difference derives from the diffusive correction 
implemented by CosmoRec vl.3 being calibrated at fixed Y p = 0.24. 

above the level of a few times 0.1%. As mentioned in the introduc- 
tion, recent constraints on the helium mass fractions from CMB 
alone derived using WMAP and ACT data dPunklev etaUl201lh 
indicate a slightly larger helium abundance of Y p = 0.313 ± 0.044. 
Varying the helium mass fraction, we obtain the modification 
shown in Fig. Q7] For all curves except the last electron scatter- 
ing is neglected. The main reason for the discrepancy is that the 
precomputed correc tion function which mimic s the effect of line 
scattering diffusion jRubino-Martm et alj2008h in the escape prob- 
ability from the 2'Pi - 1 1 So line is calibrated at fixed Y p = 0.24. 
Apparently CosmoRec vl.3 overestimates this correction slightly 
for larger values of Y p . Electron scattering only affects the recombi- 
nation dynamics at very early times, but does not have a significant 
effect on the CMB power spectra. 

We computed the modifications for the CMB power spectra 
and found that for Y p = 0.4 the difference between CosmoRec vl.3 
and CosmoRec v2.0 exceeded |AC//Q| - 0.05% only at / > 3000. 
The difference related to a full treatment of electron scattering was 
negligible. This indicates that the improvements of this work re- 
main small, and even for rather extreme values of Y p , CosmoRec 
vl.3 can be reliably used to compute the recombination history. 

The importance of line scattering alone is shown in Fig. 1181 
we internally compare the results obtained with CosmoRec v2.0 
switching the corresponding terms on and off. HyRec does not ac- 
count for the diffusion correction, however, as Fig.[T7]shows, when 
moving further away from the standard cosmology the differences 
increase to ANJN C =* -0.63% at z - 1780 for Y p = 0.4. Still 
the modifications in the CMB power spectra turn out to be only 
marginally important at high I (cf. Fig.|19t. since the impact of he- 
lium recombination on the power spectra is much smaller than the 
signals generated close to maximal visibility at z — 1 100. 

3.5.2 Direct comparison with Recfast vl.5 

Finally we compare the results for the CMB power spectra com- 
puted using the helium radiative transfer module of CosmoRec v2.0 



Effect of line scattering on the recombination history 
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Figure 18. Total correction to the ionization history for the standard cos- 
mological model with varying helium abundance. This figure illustrates the 
importance of the line scattering, as computed with CosmoRec v2.0. The re- 
combination history with line diffusion deactivated was used as reference. 

with those obtain ed using Recfast vl.5, which is the current 
standard in Camb dLewis et alj|2000h . The differences are shown 
in Fig. [20] where we again vary the helium abundance. We 
use the default parameters for the fudge functions of Recfast 
vl.5. These were ca l ibrated on the recombination calculations of 
Rubino-Martm et H] d2010h . However, even for Y p = 0.24 small 
differences in the CMB power spectra are visible at high /. These 
increase for larger values of Y p , reaching = 0.4% at / = 1000 for 
Y p a 0.4. This shows that when a wider range of parameters is 
allowed, simple fudging does not represent the detailed recombina- 
tion physics well. In this case a detailed recombination code should 
be used. Alternatively, an iterative pro cedure that relies on re cali- 
bration of the fudge factors is possible dShaw & Chluball201 ll) . 

For clarity, we emphasize that the differences seen in Fig. [20] 
are not mainly because of the improved helium recombination 
model discussed here, but were already present in comparison with 
CosmoRec vl.3. Furthermore, we would like to mention that the 
final code comparison is still not completed, so that the final cali- 
bration of the Recfast fudge factors for the current best fit cosmo- 
logical model cannot be reliably performed at this point. 

We also mention that the differences found here with Recfast 
are responsible for only a very small part of the actual difference 
seen when going from Y p = 0.24 to Y p = 0.313 (the TT power 
spectrum amplitude changes by - 10%). Therefore the recombi- 
nation corrections discussed here cannot cause the apparent slight 
tension of the ACT and SPT results for Y p with standard Big Bang 
nucleosythesis. 

3.6 Effect of electron scattering during H i recombination 

For completeness, we also computed the effect of electron scat- 
tering during hydrogen recombination. Figure [2T] presents the re- 
sults. We include all two-photon and Raman-processes up to n = 8 
in the computations. The solid black line shows the correction us- 
ing the Fokker-Planck approximat ion, which was already computed 
earlier dChluba & Sunvaev 20093). The red/dash-dotted line shows 
the result obtained using the Doppler kernel for the redistribution 
function. This kernel does not include electron recoil and Doppler 
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Figure 19. Corrections in the CMB TT (upper panel) and EE (lower panel) 
power spectra caused by neglecting the line diffusion correction. We use the 
standard cosmological model but vary the h elium abundance. T he dotted 
lines indicate the 3/1 benchmark suggested bv lSeliak et all 120031) . 



Figure 20. Corrections in the CMB TT (upper panel) and EE (lower panel) 
power spectra in comparison with Recfast vl.5. We use the standard cos- 
mological model but vary the h elium abundance. T he dotted lines indicate 
the 3 // benchmark suggested by |Seliaketai]j2003l) . 



boosting, however, as one can see these terms lead to a rather dif- 
ferent overall behaviour of the correction (blue/dashed line, which 
is computed using the Kompaneets kernel that includes those pro- 
cesses in lowest order). In particular electron recoil is important to 
get the total modification right. Still the effect of electron scattering 
can be completely neglected for the analysis of future CMB data. 

Our result is very similar to the one obtained earlier by 
lAli-Haimoud et alj feOldh . however, here we include the term due 
to Doppler boosting in addi tion, which leads to a sma ll correction. 
Furthermore, in the work o f lAli-Haimoud et alj < l20ld) electron re- 
coil is included by enforcing detailed balance rather than using the 
corresponding terms in the kernel hv/m c c 2 ). We checked that 
higher order relativistic co rrections provided by the formulae of 
ISazonov & Sunvaevl d2000h did not affect the result at any signifi- 
cant level (both for hydrogen and helium recombination). 



3.7 Effect of quadrupole lines during H i recombination 

The effect of nd-ls quadrup ole lines during hydro gen recombina- 
tion has been considered bv lGrin & Hiratal bold) . Since the tran- 
sition frequencies of the ud-ls quadrupole lines and np-ls Lyman- 
series practically coincide, it is clear that photons emitted in a nd- 1 s 
quadrupole transition are directly reabsorbed by the optically thick 
Lyman-« resonance. Similarly, photons released in a Lyman-ra tran- 
sition allow re-exciting electrons to the «d-state. 

The overall distortion introduced by nd-ls quadrupole tran- 
sition is much weaker than the one from the Lyman-n resonance. 
Therefore the photon distribution around the transition frequen- 
cies v„ )s is approximately given by the distortion introduced by the 
Lyman-n resonance. Since the nd-ls quadrupole lines are much less 
optically thick than the rcp-ls lines, the mean occupation number of 
photons in the nd-ls quadrupole lines is given by 11 !/; ss N^/SN^, 
where N„ p and Ni s denote the populations of the np and Is level. 

11 Because the variation of the photon occupation number across the nd-ls 
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Figure 21. Effect of electron scattering on the ionization history during 
hydrogen recombination. See Sect. l3.6l for details. 
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Figure 22. Effect of quadrupole lines on the ionization history during hy- 
drogen recombination. See Sect. l3.7l for details. 



Therefore the net transition rate in the rcd-ls quadrupole line can be 
approximated as AR„ ilf , = A„ ils [N ni -5Ni s n\ ~ A n6ls [N„ i -5/3N nv ], 
where A n a s is the ;id-ls transition rate, which is typically of order 
10 2 s _1 for small n (see Table lBlt . 

This is the s ame approximation that was given by 
I Grin & Hiratal j201fj|) and it effectively results in a direct nd- 
np transition. We confirm this approximation by including nd- 
ls quadrupole lines into the radiative transfer computation. We 
find that in the nd-ls series, the main effect is due to the 3d- Is 
quadrupole line, with the approximation given above being in ex- 
cellent agreement with the full numerical result. 

In addition, we consider the effect of quadrupole lines among 
excited levels. The rates can be computed as explained in Ap- 
pendix [B] We compute effective recombi nation rates (applying the 
method of I Ali-Haimoud & Hiratal [20101 to helium) including all 
quadruple lines with n < 300 (n < 100 is sufficient for conver- 
gence). We find that close to the maximum of the Thomson visibil- 



quadrupole resonance is very small the integral over the line profile can be 
carried out directly. 



ity function the quadrupole lines among excited levels are actually 
more important than the «d-ls transitions, however, the overall ef- 
fect is completely negligible (see Fig. I22t. 



4 CONCLUSIONS 

We have developed an improved helium radiative transfer module 
for CosmoRec, accounting for several subtle processes that previ- 
ously have been omitted or treated more approximately. This mod- 
ule is now included by CosmoRec v2.0, and can be activated as 
needed for most settings without affecting the runtime significantly. 

For the standard cosmology, the processes described here lead 
to small corrections in the ionization history at z — 2000, and 
the associated modifications in the CMB power spectra are neg- 
ligible (see Sect. 13.5b - However, when considering a wider range 
of cosmological parameters, the new helium recombination mod- 
ule allows more precise computation of the ionization history. For 
example, for Y p = 0.4 the modification of the ionization history 
reaches \AN e /N e \ ^ 0.22% at z - 1800, however, the overall 
correction to the CMB power spectra is small, remaining below 
|ACf/C/| - 0.05% at / < 3000. This means that the improvements 
of CosmoRec v2.0 with respect to CosmoRec vl.3 are negligible. 
We furthermore showed that the effect of line scattering during he- 
lium recombination is important only for large Y p (cf. Sect. 13.5. it . 
We also compare the results obtained with CosmoRec v2.0 to those 
computed using Recfast vl.5 and find that the parameters for the 
fudge functions probably should be re-calibrated once the final re- 
combination code comparison has been completed (see Sect. |3.5.2l . 

Finally, we considered the effect of quadrupole lines during 
hydrogen recombination (Sect. 13. 7t and found that these can be 
safely neglected in computations of the CMB power spectra. 
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APPENDIX A: RECURSIONS FOR REDUCED DIPOLE 
MATRIX ELEMENTS 

The dipole transition rates ( selection rule A/ = ±1) fo r hydrogenic 
atoms are given by (e.g., see IStorev & Hummer!ll99lh 



7 4 a*c_ 
Oflo 



,(£1) 



max(Z,/') .a) v „i ,2 

~2iTT l x »'<' 1 



(Al) 



Here Z is the atomic charge; ^ rcd = [1 + m e /M A ]~ [ is the reduced 
mass of the atom with mass Ma divided by ra e ; a is the fine struc- 
ture constant and c the speed of light; a denotes the Bohr radius 
(with out reduced mass effect) and {[) X'^, I , is the reduced matrix el- 
ement of the transition between the levels (n, I) — > (n', /'). 
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Fractional contribution of E2 lines to the lifetime of the level 




Figure Al. Fractional contribution of electric quadrapole lines to the total 
(vacuum) lifetime, 1 /T, of levels for different shells in hydrogen. Here Te2 
is the contribution to the decay rate of the level from electric quadrapole 
lines alone, while T t denotes the total decay rate. 



Table Bl. All spontaneous electric quadrapole transition rates of hydrogen 
for n < 5 and n' < n. 
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The recur sions for the reduced matri x elements of hydrogenic 
atoms given bv lStorev & Hummed J 1 99 ll> can be cast into the form 
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where we defined C\ = in 2 - P/n. These recursions can be started 
at /' with initial conditions (for n> n' > 1) 
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APPENDIX B: RECURSIONS FOR REDUCED 
QUADRUPOLE MATRIX ELEMENTS 

The electric quadrupole transition rates (selection rule Al = 0, ±2) 
for hydrogenic atoms are given by jHevl200rj : lGrin & Hiratal201fjh 
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where X"^, is the reduced quadrupole matrix element. 



Bl Case A/ = 



According to Hey ( 2006) the matrix elements <2) X' i ! ,' / are determined 
by the recursion formula. PI 



c?; i (2) ^/, + + 1 1 +2(/+i) ( "x;; 



f 



C" 



(B2) 



with the initial condition {1) X"J' n , = 0, assuming that the dipole ma- 
trix elements are known. Electric quadrupole transitions between 
s-states are forbidden. 



12 We absorb the factors ao/Z in the definition of the matrix elements 
© 0000 RAS, MNRAS 000, 000-000 



B2 Cases Al = ±2 

Following lHevl ( 12006)) the matrix elements (2) X^,_ 2 and (2) ^,' /+2 are 
determined by the recursion formulae 

(2i - i)c M a) x;Jt\ +2C( {2) X»J, 



(2) x „l 



(2) yn ; 

A n'l+2 



(21+ dc?;, 

2(2/- l)(3/ + 2)<"X;/,_, 

+ (2i+ i)q;, ' 
(2/ + 3) cf +3 <2) x;;/ / y 3 + 2 c/ +2 ^x*, >£ 

(2l + 5)q +l 
2(2/ + 3)(3/ + 8) <1) X"' + ' 



(B3) 



(2/ + 5)C" 



(B4) 



with the initial condition (2) X", n ';\ = -2n' (C'^JC,) {2) X n, /r\ , as- 
suming that the dipole and diagonal qu adrup l e matr ix elements are 
known. We note that in comparison to iHevI d200fj|) we multiplied 
all matrix elements by the phase factor (-1)"~" -1 . The Al = 
quadrupole matrix elements are in phase with this, however, for 
the off-diagonal elements a factor of (-1) remains. 

In Table lBll we give the electric quadrupole transition rates of 
hydrogen for n < 5 and n' < n. Transitions between states with 
1 = and /' = are not allowed. The typical fractional contribution 
of quadrupole transitions to the total decay rate of the initial levels 
is about 10~ 7 - 10~ 5 for excited states with n » 1 and / < 30 (see 
Fig. lAll l. The lifetimes of d-states are affected by a 10~ 5 , while s- 
and p-states are only modified at the ^ few x 10~ 7 . Beyond I a 30 
the effect of quadrupole lines on the total decay rate drops below 
10~ 7 . For highly excited levels at large / and I <n n the relative 
contribution to the total lifetime falls off like = V 1 . 
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